Temporal chaos versus spatial mixing in reaction-advection-diffusion systems 
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We develop a theory describing the transition to a spatially homogeneous regime in a mixing flow 
with a chaotic in time reaction. The transverse Lyapunov exponent governing the stability of the 
homogeneous state can be represented as a combination of Lyapunov exponents for spatial mixing 
and temporal chaos. This representation, being exact for time-independent flows and equal Peclet 
numbers of different components, is demonstrated to work accurately for time-dependent flows and 
different Peclet numbers. 
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Complex spatiotcmporal dynamics has attracted large 
interest in the last decades. Recently, reaction-diffusion 
equations have been subject of intense research due to 
their rich variety of patterns. They describe many im- 
portantphysical systems, such as chemical reactions 
lasers 0, or semiconductors On the other hand, 
complex structures can be created in fluid mechanics by 
spatial mixing In this Letter, we consider the com- 
bined action of the above effects, leading to a reaction- 
advection-diffusion system, see Eq. (0 below. Such sys- 
tems have been investigated with respect to front prop- 
agation excitable dynamics E|, and the filamental 
structure of reactive particles y|. Practically, stirred 
flows with reaction are relevant from large scales (plank- 
ton dynamics in the oceans ||) to microscales (construc- 
tion of a lab on a chip ||), and are important for bio- 
physical, ecological, and chemical applications; similar 
equations describe the dynamo effect in magnetohydro- 
dynamics |Toj . 

In this Letter, we consider a temporally chaotic reac- 
tion process. It is known that in presence of diffusion, 
temporal chaos can lead to the appearance of nontrivial 
spatial structures and space-time chaos. We demonstrate 
that such structures can appear in the presence of mixing, 
too. We develop a theory for the transition from spatially 
homogeneous (fully mixed) temporally chaotic state to a 
nonhomogeneous one, and compare it with calculations. 
Our approach is strongly related to the theory of com- 
plete synchronization of coupled chaotic systems, which 
is significantly extended because the spatial mixing leads 
to special types of coupling. 

We now formulate the reaction equations. The evolu- 
tion of the concentrations fa, i — 1, . . . , N due to reaction 
is described by a nonlinear system 



d<t>j 
dt 



% (1) 

with regular or chaotic solutions $(t). Additionally, 



each component is subject to diffusion and advection by 
an incompressible velocity field v(r, t). Normalization by 
the characteristic advection time allows a description by 
dimensionless diffusion constants dk (equivalent to Peclet 
numbers d~ l ~ Pei, generally different) and the dimen- 



sionless reaction rate, or Damkohler number, Da. 
resulting spatio-temporal equations are 



dt 



+ (v • V)fa = d^fa + DaF,(0i, 



The 



(2) 



We assume that the concentrations do not influence the 
flow, so that the field v(r, t) does not depend on fa. We 
supply Eq. (01 with no-flux boundary conditions Vfa\s = 
on the boundary S (periodic boundary conditions are 
straightforwardly treated in an analogous way). 

For a homogeneous spatial distribution of concentra- 
tions fa the advection and diffusion term vanish and the 
equations (J2J) reduce to with rescaled time. To study 
the stability of spatially homogeneous solutions 0° (Da t) 
of (0 we linearize the equations near this solution and 
obtain for a small perturbation field ipi(r,t) 



^ + (v V)cpt = diV 2 <Pi 
at 



Da Cij(D&t)<fij 



(3) 



where Cij(t) 



dFi 



is the Jacobi matrix of the sys- 
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tern Q on the solution [fa°(t), . . . , 4>%{t)]. Generally, 
the solutions of J2J grow or decay exponentially in time 
I \(p\ I ~ e At , where A belongs to the spectrum of Lyapunov 
exponents (LE). Clearly, the LEs of the solution of (JTJ 
belong to this spectrum, describing growth or decay of 
homogeneous perturbations. The stability of the homo- 
geneous solution towards inhomogeneous perturbations is 
described by the largest LE corresponding to a spatially 
varying Lyapunov vector ip(r, t) - the transverse LE Aj_ 
[in the sense that it is transverse to a manifold of spatially 
homogeneous solutions of QJ]. For diffusion constants 
and time-dependent flow v(r, t) given, the transverse LE 
can be determined only numerically. For certain situa- 
tions, considered below, we obtain this exponent analyt- 
ically. 

We start the analysis of ® with the simplest case, 
a time-independent velocity field v = v(r) and equal 



2 



diffusion constants di — d. In this case the ansatz 
ipi(r,t) = X(r)$>i(t) allows for a separation of time and 
space dependence of the perturbation field. The spa- 
tial component is determined by the advection-diffusion 
eigenvalue problem 



dV 2 X - (v • V)X = -TX, 



VXI 







(4) 



The eigenvalue T describes the decay of non- 
homogeneous states of the passive scalar field u(r, t), gov- 
erned by 

^ + (v • V)u = dV 2 u , (5) 

because for an exponentially time-dependent solution, (J5J 
reduces to (@J. This problem has been recently analyzed 
in The eigenvalue T = corresponds to the spa- 

tially homogeneous solution and does not contribute to 
the stability of inhomogeneous perturbations. For the 
latter, the mode corresponding to the smallest non-zero 
eigenvalue of Q is relevant, we denote it by 7. As has 
been argued in [llj], this eigenvalue crucially depends on 
the nature of flow v(r), and thus on Pe. If the flow 
is mixing, the eigenvalue 7 is well separated from zero 
even for small diffusion constant d. However, a flow typ- 
ically contains chaotic and regular domains (islands of 
Kolmogorov-Arnold-Moser tori). Then, for small diffu- 
sion constants, there are weakly decaying modes concen- 
trated in these islands, so that 7 may be rather small. We 
notice, that in terms of effective diffusion (coarse grained 
on the system size L) this eigenvalue can be represented 
as 7 = d e ffL~ 2 [l2j| . 

The equation for the temporal part $i is 



•DaCy(Dai)$ 3 - 



(6) 



With the ansatz <I>i = e~~ ' Wi this equation is trans- 
formed to the equation for linear perturbations of the 
reaction problem (JIJ 

— - — Da, Cij(Da,t)wj , 



dt 



(7) 



its asymptotic solution is ||w|| ~ e *, where A is the 
largest LE of the attractor in QJ. Thus, for the pertur- 
bation we have <I>i ~ e ( DaA -7)* anc i the explicit formula 
for the transverse LE reads: 



A 1 



Da • A - 7(Pe) 



(8) 



The stability condition of the spatially homogeneous 
state can be formulated as Aj_ < 0. If the oscillations 
of the concentrations are regular, then the largest Lya- 
punov exponent is nonpositive A < and this regime is 
always stable against spatially inhomogeneous perturba- 
tions. A nontrivial transition occurs for a chaotic regime, 
if A > 0. Here the stability condition leads to the critical 
value 

7 (Pe) 



Da c 



A 



(9) 



A similar condition for a trivial case of reaction-diffusion 
system has been obtained in Il3l and for an abstract 
mapping model of mixing in |l4| . We emphasize that 
condition 0, obtained for a realistic reaction-advection- 
diffusion system, can be directly applied to an experi- 
ment. Indeed, the eigenvalue 7 can be directly measured 
from the time evolution of the contrast of a passive scalar 
in the flow under investigation 15]. The LE A can be 
determined from the advection-free setup: the critical 
domain size L c at which the patterns appear is related 
to A via CdL~ 2 — A, where C is a geometrical factor 
depending on the domain form. 

The analysis above is based on simplifying assump- 
tions: time-independence of the velocity field and equal- 
ity of diffusion constants. In general, if the velocities are 
time-dependent and the diffusion constants are different, 
Eq. cannot be simplified and should be analyzed nu- 
merically. Since we are interested in stability with re- 
spect to spatially inhomogeneous perturbations, the so- 
lution should be sought in the class of fields tp having zero 
spatial average (the spatially homogeneous modes in 
are decoupled from other modes). Numerically, one can 
use the usual method for calculation of the largest LE: 
starting with an arbitrary initial field in (0, with van- 
ishing spatial average, one integrates 0) along with 
performing normalization of the linear field to avoid nu- 
merical over- or underflow; averaging the logarithm of 
the normalization factors yields the transversal LE Aj_. 

We apply this numerical method to the time-dependent 
flow, 2-7r-periodic in space, suggested in (lrj| : 

v = K/(t) cos[y + e x (t)], Oy[l - f(t)] cos[x + 9 v (t)}} . 

(10) 

Here, the simple analytical analysis above is not appli- 
cable. Depending on the functions f(t),9 x (t),9 y (t), the 
flow can be time periodic or irregular. In the former 
case, the particle trajectories demonstrate typical Hamil- 
tonian dynamics with chaotic regions and stability is- 
lands (see 0, ITjj)- A weakly turbulent irregular flow 
occurs if the functions f(t), 9 x (t), 9 y (t) are random func- 
tions of time. We have considered both of these cases 
with the reaction dynamics given by the Lorenz equa- 
tions. 

An irregular flow was modeled by setting f(t) as a 
(0, l)-telegraph process with independent exponentially 
distributed time intervals Ti nt , and independent uni- 
formly distributed phases 9 X , 9 y £ (0, 2w). The transverse 
Lyapunov exponent has been calculated as described 
above; the results are shown in Fig. ^ Remarkably, Aj_ 
is nearly a linear function of Da like in 0. Moreover, 
we can demonstrate numerically that the formula is 
valid also quantitatively when all the diffusion constants 
are equal. To this end, we have calculated the decay rate 
7 from the linear evolution Eq. of the passive scalar. 
Now this equation does not yield the simple eigenvalue 
problem Q because the velocity field is time dependent, 
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considers two coupled nonlinear systems 



FIG. 1: Numerical calculation of transverse LEs from Eqs. 
© and from approximations JHJ and 1161 1. Parameters: a x = 
a y = 20, (Tint) = 0.5. Squares: di = 0.1; circles: di = 0.1, 
d,2 = 0.2, ds = 0.5. Lines are according to (0 and l|16|l . 



but defines in full analogy with the discussion above the 
asymptotic decay rate 7 in the sense of the LE (cf. 0]): 



7 



= — lim 



In I Kr,t) 1 1 



t 



(11) 



Thus, 7 is physically interpreted as the asymptotic decay 
rate of the contrast of the passive scalar in the advection- 
diffusion problem In Fig. ^ we compare the trans- 
verse LE A^, calculated from JBJ), with the numerical 
estimation of A^. Here, the value of 7 has been calcu- 
lated according to i|ll|) , A is the largest LE of the Lorenz 
model; the correspondence with the numerics is within 
statistical errors. This result indicates that the time de- 
pendencies due to chaotic time evolution of species and 
due to irregularity of the flow are essentially "separable." 

The analysis above must be modified if the diffusion 
constants in © are different. In this case one has differ- 
ent passive scalar evolution problems Eq. © , and Eq. © 
is modified to 



dt 



-7i$j + DaCy(Dat)^ . (12) 
The generalized transverse LE for (|12|) is 

A_i_(Da,7i, ...,7iv) , (13) 

for 71 = . . . = 7tv = 7 it reduces to but generally 
cannot be related to the LE of the reaction system like 
in J7|). As a result, instead of © and 10, we obtain the 
following condition for the critical Damkohlcr number 



Ax(Da cr ,7i, . . . ,7jv) = . 



(14) 



There is a connection of the defined above transverse 
LEs with the theory of synchronization. In the latter one 



dt 
dt 



(15) 



and looks for the stability of the completely synchronized 



state <f>t(t) = 4>Y'{i)- Then the linearized equation for 
perturbations coincides with (|12fl . Thus, the generalized 
transverse LE Aj_ determines the synchronization thresh- 
old. This analogy shows that in a stirred reaction the role 
of coupling constants 7, is played by effective decay rates 
where both advection and diffusion contribute. 

In general, some coupling constants can be absent, i.e. 
the corresponding coefficients 7, vanish. In the context 
of chaotic mixing, such a situation appears, e.g., for sur- 
face reactions. Here, it can happen that mobility of some 
chemical species is large so that they are advected by the 
fluid flow, while other species are so strongly chemisorb 
that they cannot move laterally across the surface (see, 
e.g., 01). These chemisorbed species are described by 
vanishing spatial decay In the biological context, 
mixed and nonmixed components interact, e.g., in ma- 
rine sediments by tidal flows. 

Notice, that although a transformation <I>i = e~ st Wi 
does not allow us to reduce (|T2*)l to J7J, one parameter, 
sa y In, can be eliminated with such a transform. As a 
result, the generalized transverse LE l|13l) can be repre- 
sented as a function of N — 1 parameters only: 



Aj_(Da,7i) = DaAj_ 1 



7i ~ In 
Da : 



In-i - 7n 
Da 



, 0j -7jv 
(16) 

In Fig. 13 we show A^ for the Lorenz system. The above 
result for A^ agrees well with the direct calculation of 
transverse LEs for the full system , see Fig. ^ 




FIG. 2: Generalized transverse LEs for the Lorenz model ac- 
cording to ll 1 12 tt with 7z = and Da = 1. The zero contour 
level is shown with a line. 

Above we have performed the stability analysis of the 
spatially homogeneous regime. Next we discuss prop- 
erties of structures that appear beyond the stability 
threshold. A natural quantity to characterize the in- 
homogeneity of the pattern is the contrast, or variance, 
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W. t = <(0i-(0i)) 2 ). We solved © using fT® with the pe- 
riodic function f(t) = (0, 1) like in [16( with unit period 
and constant 6 x ,yi an d the Lorenz model as reaction. All 
three components show the same onset, consistent with 
linear theory, see Fig. |3| Near the onset of spatial inho- 
mogeneity, the temporal behavior of the pattern contrast 
is highly intermittent (see inset in Fig.|3J). As in a transi- 
tion to complete synchronization, the reason for this in- 
termittency is the fluctuations of the local exponents Aj_ 
and 7^, similar to previously investigated cases 0>E1- 




Concluding, we have applied the method of general- 
ized transverse Lyapunov exponents to the analysis of 
transition from homogeneous to inhomogeneous field in 
a stirred, temporally chaotic chemical reaction. Our the- 
ory is valid for general chaotic reactions and general mix- 
ing setups - the latter do not need to be perfect (hyper- 
bolic). For equal diffusion constants, the critical value of 
the Damkohler number © has a clear physical meaning: 
it determines, whether the growth rate due to temporal 
chaos (given by the maximal LE of the nonlinear sys- 
tem) wins over the decay rate due to mixing (given by 
the decay rate of the contrast of a passive scalar in the 
flow). For different diffusion constants, or for surface re- 
actions where only some species are stirred, the critical 
Damkohler number can be formulated as a novel problem 
in the theory of complete synchronization, whose solution 
is expressed in terms of generalized transverse LEs. 

We acknowledge fruitful discussions with B. Eckhardt, 
O. Popovych, T. Tel, and G. Zaslavsky. A. S. thanks 
DAAD for support. A. S. and M. A. were supported by 
the German Science Foundation. 



FIG. 3: Time- averaged contrast of the reaction species vs the 
Damkohler number [the critical value is is 1.096, cf. (j5)]. The 
inset shows temporal intermittency of the contrast W\ for 
Da = 1.15. Parameters: di = 0.04 • it 2 , a x = a y — 1.2 ■ it. 

We have demonstrated that the transition to spatially 
inhomogeneous structures in a mixed flow with chaotic 
in time evolution of concentrations is determined by the 
transverse LE. This exponent can be with good accuracy 
represented by a sum of "temporal" and "spatial" contri- 
butions according to ©, or, more generally, to (jTT)|l . We 
can give the following physical arguments in favor of this 
effective "decoupling" allowing for the separation ansatz. 
(i) For large Peclet numbers, although different species 
have different diffusion constants, their mixing properties 
are to a large extent determined by advection. Thus, one 
can expect that spatial structures of Lyapunov vectors for 
different variables ipi(r,t) in Q are close to each other. 
This is confirmed by numerical simulations: we have cal- 
culated the spatial correlation coefficients between differ- 
ent variables for the example presented in Fig. ^ (the case 
of different diffusion constants) and have found for these 
correlations the values larger than 0.83. The same prop- 
erty can be expected for small Peclet numbers, where the 
spatial structure is given by the purely diffusive version of 
and one has the same modes for different species, only 
the decay constants are different, (ii) Although stirring 
for a time-dependent flow is not constant in time, this 
seems to have a very small effect on the average trans- 
verse LE. We have checked this by solving the generalized 
system Ijl2|l where the decay rates 7, are not constants 
but oscillating functions of time. In some range of peri- 
ods, the corresponding LE Aj_ remained the same within 
1% even if the modulation was as large as 100%. 
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